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ABSTRACT 

A  method  for  identifying  and  -  after  identification  - 
reducing  linear  time  invariant  systems  on  the  basis  of  input  - 
output  records  of  the  system  is  reviewed,  applied  to  physical 
devices  and  extended  to  handle  the  case  where  noise  is  present 
The  method  is  implemented  using  an  analog  to  digital  converter 
in  order  to  have  the  input  and  output  records  in  a  digitized 
form  and  thus  to  be  able  to  use  a  digital  computer  program 
composed  of  a  numerical  integration  subroutine  and  another 
subroutine  to  solve  overdetermined  sets  of  linear  algebraic 
equations.   Several  practical  examples  are  presented  to 
demonstrate  the  accuracy  and  present  capabilities  of  the 
procedure. 
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I.   INTRODUCTION 

A.   THE  IDENTIFICATION  PROBLEM 

All  the  modern  techniques  of  analysis  and  design  of 
any  kind  of  system  are  based  on  the  previous  knowledge  of 
a  mathematical  model  of  the  system  to  be  studied.   This 
principle  is  also  applicable  to  control  systems.   The  form 
of  the  model  depends  on  the  methods  of  analysis  and  design 
to  be  employed  as  well  as  on  the  physical  characteristics 
of  the  system  and  since  most  of  the  theory  on  control 
systems  is  based  either  on  the  state  space  or  transform 
representation  of  systems,  the  majority  of  the  mathematical 
models  will  be  a  set  of  state  equations  or  a  transfer 
function. 

After  deciding  what  basic  form  the  mathematical  model 
should  take,  the  problem  of  determining  the  numerical  values 
of  the  parameters  arises.   These  values  can  often  be  deter- 
mined from  the  laws  of  physics  and  the  data  supplied  by  the 
manufacturer  or  obtained  through  testing.   This  is  not  always 
the  case,  however;  sometimes  the  laws  of  physics  become 
mathematically  intractable  or  are  very  difficult  to  apply  and 
quite  often  the  values  of  certain  key  parameters  are  not 
available  or  have  changed  with  time.   It  is  in  these  cases 
when  the  identification  problem  arises.   This  problem  consists 
in  determining  the  system  model  given  the  input,  output  and 
initial  conditions. 


For  the  purpose  of  this  work,  the  identification  problem 
can  be  stated  as  follows: 

Given  a  system  considered  a  black  box,  obtain  a  record 
of  the  input  and  output  of  the  system  over  some  finite 
period  of  time  -  400  msecs.  in  this  particular  case  -  and 
from  it  find  a  mathematical  model  -  the  transfer  function 
relating  output  and  input  in  this  case  -  and  the  values  of 
the  parameters  in  such  a  way  that  the  model  will  accurately 
describe  the  system. 

It  should  be  kept  in  mind,  however,  that  the  problem  of 
identifying  a  system  solely  on  the  basis  of  input  and  output 
quantities  is  very  rare  in  engineering,  because  in  most  of  the 
cases  the  engineer  has  some  idea  whether  the  system  is  linear 
or  non-linear,  time  varying  or  time  invariant,  what  is  the 
order  of  the  system  itself  and  perhaps  what  are  the  time 
constants  involved  and  in  the  presence  of  noise  the  statis- 
tical characteristics  of  it.   For  these  reasons  rather  than 
having  a  BLACK  BOX,  what  the  engineers  have  is  what  could  be 
called  -  a  GREY  BOX  -,  and  as  will  be  shown  later,  the 
identification  techniques  depend  very  heavily  on  knowledge 
of  the  characteristics  mentioned  above  and  in  this  particular 
case  THE  ORDER  OF  THE  SYSTEM. 

Although  the  control  systems  literature  on  system  identi- 
fication techniques  is  quite  vast,  there  are  no  known  identi- 
fication schemes  capable  of  handling  all  identification 
problems  effectively.   Some  of  the  proposed  identification 
methods  are  mentioned  in  the  following: 


Nieman,  Fisher  and  Seborg  [1]  have  pusblished  a  paper 
summarising  most  of  the  common  approaches  to  the  problem  of 
identifying  lumped  parameter  systems.   Eykhoff  [2]  has  a  good 
discussion  of  the  industrial  application  of  several  methods. 
Since  the  response  to  an  impulse  gives  directly  the  transfer 
function  of  linear  systems,  obtaining  this  response  is  a 
common  scheme  used  in  linear  identification.   Mishkin  and 
Haddad  [3]  have  developed  a  technique  for  finding  the  impulse 
response  based  on  samples  of  the  system  output  and  its  deri- 
vatives.  Levin  [4]  and  Kerr  and  Surber  [5]  have  developed  a 
scheme  for  estimation  of  the  impulse  response  from  the 
samples  of  a  noisy  input  and  output.   Turin  [6]  and  Lichtenberger 
[7,  8]  have  used  a  matched  filter  to  obtain  an  identification. 
Hill  and  McMurtry  [9]  have  suggested  the  use  of  white  noise, 
a  binary  test  signal  and  crosscorrelation   techniques  to 
attack  the  problem.   Cooper  and  Lindenlaub  [10]  have  investi- 
gated the  use  of  matched  filter  and  crosscorrelation. 

Another  common  approach  to  the  problem  is  to  determine 
the  coefficients  of  the  differential  or  difference  equation 
which  describes  the  system.   Kumar  and  Sridhar  [11]  have 
employed  the  method  of  quasilinearization  with  some  success. 
Nagumo  and  Noda  [12]  have  developed  a  learning  approach  to 
the  problem.   Bass  [13]  uses  modulating  functions  and  the 
method  works  well  in  the  presence  of  noise.   Astron  and 
Bohlin  [14]  use  a  statistical  optimal  method  of  determining 
differential  equation  parameters  known  as  the  "maximum 
likelihood  method".   Peterka  and  Smuk  [15,  16]  have  developed 


a  similar  method  which  is  not  optimal  but  is  considerably 
simpler  computationally  and  it  is  known  as  the  "prior 
knowledge  fitting  method".   Ho  and  Kalman  [17]  have  proposed 
an  algorithm  for  determining  state  variable  models  of  sampled 
data  systems  which  works  well  in  the  presence  of  noise,  and 
has  been  extended  to  continuous  systems  by  Eldem  [18] . 

Identification  methods  vary  widely  with  respect  to  how 
much  must  be  known  about  the  system  before  the  method  can  be 
applied.   In  general,  the  less  known  about  a  given  system, 
the  more  difficult  it  is  to  find  a  method  capable  of  accom- 
plishing the  identification. 

B.   OBJECTIVES  OF  THIS  PAPER 

This  paper  will  present  a  research,  handling  actual 
hardware,  of  an  identification  technique  originally  suggested 
by  Diamesis  [19]  and  studied  using  simulation  techniques  by 
W.  R.  Hansell  [20] ;  also  an  attempt  at  using  the  technique 
in  the  presence  of  noise  is  made.   The  method  requires  a 
knowledge  of  the  system  input  and  output  over  some  finite 
time  interval  and  a  rough  estimate  of  the  order  of  the  system 
is  desired.   Although,  theoretically  the  system  input  used 
need  not  be  restricted  to  a  class  of  testing  functions,  that 
is,  it  can  be  completely  arbitrary,  experimental  results 
showed  that  inputs  such  as  step,  ramp  or  in  general  any 
periodic  wave,  sine,  square,  etc.,  should  be  avoided. 

Unlike  some  identification  schemes  which  require  the 
calculation  of  derivatives  of  the  input  and  output,  the 


technique  to  be  presented  requires  only  integrals  of  the 
input  and  output,  and  the  advantages  of  numerical  integration 
over  differentiation  are  well  known. 

In  Chapter  II,  the  theoretical  development  of  the 
identification  technique  is  given.   A  method  for  identifying 
the  initial  conditions  of  the  unknown  system  is  also  presented, 
but  in  this  case  it  is  necessary  to  know   the  value  of  the 
derivatives  of  the  input  at  time  zero.   Since  the  initial 
conditions  do  not  modify  in  any  sense  the  mathematical  model, 
their  knowledge  is  not  important. 

Chapter  III  presents  the  method  used  to  check  the 
developed  theory  in  Chapter  II.   Several  examples  are  given 
to  demonstrate  the  capabilities  of  the  method. 

In  Chapter  IV,  recommendation  and  conclusions  are  given 
and  following  this  chapter  a  complete  listing  of  all  computer 
programs  are  presented. 


II.   IDENTIFICATION  BY  MULTIPLE  INTEGRATIONS 

A.   GENERAL  APPROACH 

The  development  that  follows  is  similar  to  the  ones 
given  by  Diamesis  [19]  and  W.  R.  Hansell  [20] .   The  develop- 
ment given  by  Diamesis  is  restricted  to  the  case  where  all 
initial  conditions  are  zero  and  the  parameters  of  the  model 
are  the  unknowns  of  a  uniquely  determined  set  of  linear 
algebraic  equations.   Hansell  extended  the  development  to 
the  case  where  the  initial  conditions  need  not  be  zero  and 
he  makes  use  of  an  overdetermined  set  of  linear  algebriac 
equations  with  the  model  parameters  and  some  coefficients 
which  have  information  about  the  initial  conditions,  as 
unknowns. 

The  overdetermined  set  of  equations  are  solved  using 
the  method  of  least  squares.   The  method  is  extended  to 
handle  noisy  inputs  and  outputs. 

Any  single  input,  single  output,  lumped  parameter,  linear 
time  invariant  system  can  be  described  by  a  linear  ordinary 
differential  equation  with  constant  coefficients.    The 
general  form  of  this  equation  is  given  in  equation  (1) .   A 
set  of  initial  conditions  is  also  necessary  to  completely 


define  the  system. 

yn)  (t)  +  a^y11"1*  (t)  +  +aQy(t)  = 

(1) 
=  b  U   (t)  +  +  b  U(t)  where,  in  a  general  case 

n>m  and 
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y   (t)  and  lP  (t)  represent  the  i    and  j    derivatives 
of  the  output  and  input,  respectively. 

The  necessary  initial  conditions  are: 

The  value  at  time  zero,  of  the  output  and  its  n-1  first 
derivatives  and  the  value  at  time  zero,  of  the  input  and  its 
m-1  first  derivatives. 

The  notation  to  be  used  is: 

y(o),  y   (o)  y     (o)  ,  y     (o)  to  denote 

the,  n-1,  I.C.  of  the  output  and 

u  (o)  ,  u   (o)  um   (o)  ,  u    (o)  to  denote  the , 

m-1,  I.C.  of  the  input. 

The  identification  problem  to  be  treated  here  consists  of 
determining  the  coefficients  a.  and  b.  of  equation  (1),  on 
the  basis  of  input  and  output  records  taken  over  some  arbi- 
trary time  interval,  for  every  pair  of  values  chosen  for  n 
and  m.   The  initial  conditions  will  be  assumed  to  be  unknown 
and  also  it  is  assumed  that  the  noise  present  is  zero  mean 
and  that  the  standard  deviation  is  small  compared  with  the 
R.M.S.  of  the  input  and  output,  for  the  technique  to  work 
properly. 

Taking  the  Laplace  transform  of  (1),  yields  equation  (2) 

snY(s)  +  a   ,sn~1Y(s)  +  +  a  Y(s)  = 

n-±  o 

bmsmU(s)  +  +  boU(s)  +  Vl5""1  +  •••  +  do  (2) 

where  the  d,  coefficients  account  for  the  contribution  of 
the  initial  conditions. 

Dividing  equation  (2)  by  s    results  in  equation  (3) 
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IlSl  +a  ,  Stt2--+  +a  *I»>-b    gjgl   + 

n-1    2  o   n+1    m  n-m+1 

s  s  s         s 

...  +  b  £l|>  +d   ,-i  + +d     1 

o  sn+l    n-l^  o  _        (3) 

'Taking  now  the  inverse  transform  of  equation  (3) ,  gives 
equation  (4) . 
"t 


,.  ,  , .  n-m+1 
u(t)dt      + 


u(t)dtn+1  + 


...  (n+1)  ...!   dtn+1 

/ 

o    Jo  Jo  Jo 

And  since  the  system  is  assumed  to  be  time  invariant, 
nothing  has  been  lost  by  setting  the  lower  limit  on  the 
integrals  equal  to  zero. 

Rearranging  terms  in  equation  (4) ,  equation  (5)  is 
obtained. 
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(5) 


Vi       /  /  dt     =  -  /  *<*>« 

Since  records  of  the  input  and  output  are  assumed  to  be 
known,  a  linear  algebraic  equation  with  the  system  parameters 
a.,  b.,  and  the  terms  that  contain  information  about  the 
initial  conditions  as  unknowns,  can  be  formulated  by  performing 
the  multiple  integrations,  from  zero  to  some  time  tf,  or  from 
any  time  t   to  a  final  time  tf,  after  assuming  some  values 
for  n,  and  m.   A  set  of  2n+m+l  equations  can  be  obtained  by 
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letting  t,  take  on  2n+m+l  different  values.   Assuming  that 

the  equations  are  linearly  independent,  it  will  now  be 

possible  to  solve  for  the  2n+m+l  unknowns;  n+m+ 1 ,  a.,  b. 

coefficients  and  the  n  d.  terms  which  contain  information 

—  k 

about  the  initial  conditions.   It  will  be  shown  later  how 
those  coefficients  can  be  used  to  determine  the  initial 
conditions. 

Consider  now  the  case  when  noise  is  present  at  the  input, 
at  the  output,  or  at  both. 

Under  these  conditions  equation  (5)  becomes: 

/       ...      ry(t)+v(t)'|  dtn+1  + 


o 


fy(t)+v(t)l 


jy(t)+v(t)j 


y(t)+v(t)|   dt2  - 


bQ    /       ...  rv   ...     ;      i  ltn4 


I  u(t)+r(t)J 


|u(t)+r(t)  J 


n-m+1)  ...         |u(t)+r(t)|  dt 

(6) 


n-m+1 


d     J      ...  n+1  ...     f      dt 


n+1 
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dn_1   I         /      dt*  -  -  /       y(t)+v(t)    dt 

o       Jo  Jo 

where  v(t)  and  r(t)  represent  the  noise  at  the  output  and 
input  respectively. 

If  the  coefficients  in  equation  (6),  that  is,  the 
multiple  integrals  of  the  input,  output,  and  time,  are  going 
to  have  the  same  values  as  the  corresponding  ones  in  equation 
(5),  in  order  to  have  the  same  solution  for  the  a.,  b.,  and 
d,  ,  in  other  words,  the  same  identification,  the  presence  of 
noise  should  not  affect  the  value  of  the  integrals. 

Consider  a  general  term: 

A 

...  (p)  ...  I    ry(t)+v(t)|  dtp  (7) 


[y(t)+v(t)j 


o  Jo 

Equation  (7)  can  be  written: 

...  (p)  ...    /      y(t)dtp  +    I    ...  (p)../   v(t)dtp 


o  jo  Jo  Jo 

The  first  integral  of  equation  (8).  (8) 

...  (p)  ...   /      y(t)dtp  (9) 


is  the  true  value  of  the  coefficient  of  a    . ,  in  equation 

n-p+1     n 

(5)  and 
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(p)  .  .  . 


v(t)dtp 


(10) 


o  #  o 

the  undesired  contribution  of  the  noise.  Thus,  what  has  to 
be  done  is  to  find  when  the  value  given  by  equation  (10)  is 
negligible  with  respect  to  the  value  of  equation  (9) . 

Assuming  zero  mean  noise,  the  expected  value  of 
equation  (10)  is: 


.  . .  (p)  .  .  . 


v(t)dtp 


. .  .  (p)  .  .  . 


H 


dtp  =   0 


(id 


Since  zero  mean  noise  was  assumed  and  the  expected  value 
of  the  integral  is  equal  to  the  integral  of  the  mean  value. 
But  the  fact  that  the  mean  value  of  equation  (10)  is  zero  is 
not  enough  to  say  that  the  influence  of  the  value  of  equation 
(10)  in  (7)  is  negligible.   Besides  the  mean  of  equation  (10), 
its  standard  deviation  or  R.M.S.  value  should  be  found  in 
order  to  be  able  to  compare  it  with  equation  (9),  and  if  the 
R.M.S.  value  of  equation  (7)  is  very  small  compared  with  the 
value  of  equation  (9) ,  it  can  be  said  that  the  influence  of 
the  noise  will  not  significantly  affect  the  identification. 

The  variance  of  equation  (10)  considering,  p  =  2,  for 
simplicity,  is  given  by  equation  (12) 
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Oz   =   E 


v(Y)dydC 


(12) 


o      #  o 

where  y  and  E,   are  dummy  variables.   Equation  (12)  can  be 
written 

v(Yl)dYl  dCjL  '  I         7   v(Y2)dy2dC2 


)  v(y2)|    dYl  dY2  aZ±   d^2   (13) 


but  E  I  v  (y-j)  v(y2)  I  is  recognized  to  be  the  auto-correlation  and 

can  be  written  as  R  (t)  .  where  t  =y-.~Yt  • 

v    '  '21 

The  R.M.S.  value  of  equation  (10)  can  now  be  expressed 
as  follows: 


1 
7 


Rv(T)dY1dY2d?1dC2      < 


V(o)dY1dY2dC1dC2      \2 


^Rv(o)dY]L   dq 


rRv(o)dY2dC2 


^Rv(o)dYd? 


R    (x)<    R    (o) 
v        —     v 


(13a) 
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or  going  back  to  the  initial  notation 

ft 


^Rv(o)cH:  (14) 


fo      JO 

If  now  (14)  is  compared  with  (9),  with  p  =  2,  it  is  observed 
that  both  integrals  have  the  same  form.   Therefore,  if 

y(t)  >>  R  (o)  =  R.M.S.  value  of  v(t),  the  noise  is  not 
going  to  have  much  influence  in  the  identification.   The 
same  conclusion  can  be  reached  for  any  value  of  p. 

As  was  said  before,  the  set  of  2n+m+l  linearly  indepen- 
dent algebraic  equations  obtained  from  equation  (5)  letting 
t,  take  2n+m+l  different  values,  would  solve  the  problem  of 
identifying  the  parameters  a.,  b.,  and  d,  ,  but  since  for 
some  kind  of  inputs  and  some  t_,  the  2n+m+l  equations  may 
not  be  linearly  independent,  thus  instead  of  forming  just 
2n+m+l  equations,  a  set  is  formed  which  has  a  number  of 
equations  substantially  bigger  than  2n+m+l  and  this  over- 
determined  set  is  solved  using  the  method  of  least  squares. 

This  technique  also  allows  some  flexibility  in  the 
selection  of  the  initial  values  of  n  and  m  as  will  be  shown 
later. 

If  in  equation  (2)  the  initial  conditions  are  set  equal 
to  zero,  after  a  simple  manipulation,  the  transfer  function 
of  the  system  can  be  obtained  and  is  given  in  equation  (15). 

Y(s)     b  sm  +  .  .  .  +  b,s  +  b 

=  _m 1     o (15) 

U(s)     s   +  a   ,s    +  ...  +  a.s  +  a 
n-1  1     o 

18 


It  should  be  noticed  that  the  coefficients  a.  and  b. 
of  equation  (10)  are  exactly  the  same  as  the  a.  and  b.  of 
equation  (5).   Therefore,  as  soon  as  the  set  of  overdetermined 
algebraic  equations  is  solved  for  a.  and  b.,  the  transfer 
function  is  determined. 

Equation  (15)  is  usually  written  in  the  form  (16) . 


Y(s)    K£sm  +  cm-1sm  +  ...  +  c1s  +  cQ"j 

1,  /  \     n        m— 1 

U(s)    s   +  a   ,s    +  ...  +  a,s  +  a 
n-1  ±     o 


(16) 


Where  K  =  b   is  the  so  called  constant  gain  and  c^  =  _g_  for 


^       E3- 


m 
o  <_   g  <_  m. 

If  the  roots  of  the  numerator  and  denominator  polynomials 

are  found,  the  more  familiar  form  of  the  transfer  function 

can  be  written;  equation  (17) . 


m 


T^i    (s  +  qj} 

uT^T       5  ~   (s  +  Pi) 


(17) 


i  =  1 

Where  -p.  and  -q . ,  are  the  poles  and  zeros  respectively. 

Equation  (16)  can  also  very  easily  be  written  in  state 
form  as 


• 

0          J 

xl 

• 

X2 

— 
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y   =   I  c  c,   ....  1    0   ....  0 


xl 


x2 
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u  n 


(17a) 


But,  how  are  n  and  m  determined?  It  was  said  that  in 
order  to  apply  the  method  efficiently  a  rough  idea  of  the 
order  of  the  system  should  be  available,  although  it  is  not 
absolutely  necessary;  that  rough  idea  is  used  to  pick  n'  and 
m'  as  the  first  approximation  for  n  and  m,  taking  into  account 

that  nf>n  and  m'>m  to  be  sure  that  neither  poles  nor  zeros 
are  missing.  Since  m  is  usually  smaller  than  n,  m'  should 
not  be  greater  than  n'-l. 

Consider  then  that  n  and  m  are  the  true  numbers  of  poles 
and  zeros  and  that  n'  and  m1  are  assumed.   Under  these 
conditions  one  would  expect  to  obtain  n'-n,  or  m'-m  -  what- 
ever number  is  smaller  -  poles  equal  to  zeros,  what  is 
equivalent  to  saying  that  they  are  not  present  in  the  system, 
what  is  actually  the  case;  the  (n1  -  n)  -  (m1  -  m)  poles  or 
(m1-  n)  -  (n'  -  n)  zeros  in  excess  which  have  not  been 
cancelled  will  appear  to  be  either  in  the  right  half  plane 
or  with  values  too  big  compared  with  the  actual  values  of  the 
poles  and  zeros  in  the  system.   As  a  result  after  a  first 
guessing  for  n  and  m,  a  clear  indication  of  the  actual  order 
of  the  system  and  of  the  number  of  zeros,  should  show  up. 
Examples  in  Chapter  III  will  demonstrate  this  point. 
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B.   IDENTIFICATION  OF  THE  INITIAL  CONDITIONS 

Although  the  knowledge  of  the  initial  conditions  is 

not  important  by  itself,  it  might  be  if  one  desired  to 

compare  the  system  model  with  the  actual  system  by  exciting 

the  system  model  with  the  same  input  that  was  used  in  the 

identification.   By  comparing  the  output  of  the  model  with 

the  output  of  the  system  a  rough  idea  of  the  accuracy  of 

the  model  can  be  obtained.   However,  in  order  for  this  to  be 

possible,  all  the  initial  conditions  of  the  input  and  output 

should  be  known. 

When  the  overdetermined  set  of  linear  algebraic  equations 

is  solved,  besides  the  values  of  a.  and  b.,  the  d,  initial 
'  1      j       k 

condition  terms  are  also  found.   These  d,  terms  are  related 

k 

with  the  initial  conditions  by  equation  (18).  This  equation 
can  easily  be  found  by  specifically  writing  in  the  contribu- 
tion of  the  individual  initial  conditions  when  the  Leplace 
Transform  of  equation  (1)  was  taken,  to  obtain  equation  (2) . 
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Equation  (18)  requires  the  knowledge  of  m-1  derivatives 
of  the  input  and  numerical  techniques  may  be  necessary  to 
find  them,  however,  this  difficulty  may  be  avoided  in  many 
cases  by  choosing  the  beginning  of  the  input-output  record 
at  a  moment  where  the  input  is  constant  and  thus  all  the 
derivatives  would  be  zero. 

C.   SPECIAL  CASES  WHERE  SIMPLIFICATION  COULD  BE  POSSIBLE 

Theoretically  the  method  should  identify  the  system 
independently  of  the  number  of  unknowns, but  practically 
this  is  not  the  case  due  to  the  natural  limitations  in  the 
accuracy  of  the  numerical  integration  methods  and  in  the 
measure  of  the  input  and  output  records,  thus  it  would  be 
desirable  to  reduce  the  number  of  unknowns  as  much  as 
possible,  because  as  a  general  rule  the  fewer  the  unknowns, 
the  more  accurate  its  identification. 

A  great  reduction  in  the  number  of  unknowns  occurs  if 
the  record  of  the  input  and  output  begins  when  all  the 
initial  conditions  are  zero.   Under  these  conditions,  all 


(18) 
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the  terms  d,  of  equation  (5)  are  zero  and  the  number  of 
unknowns  is  reduced  from  2n+m+l  to  n+m+1  since  it  was 
assumed  n>m. 

Any  additional  information  could  be  used  to  simplify 
the  problem,  for  instance  if  the  value  of  the  steady  state 
gain  constant  were  known  the  number  of  unknowns  could  be 
reduced  by  one. 
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III.   IMPLEMENTATION 

A.   NUMERICAL  METHODS 

As  was  mentioned  before,  after  having  the  record  of  the 

input  and  output  of  the  system,  which  can  be  obtained  as 

described  in  Appendix  I,  it  is  necessary  to  perform  the 

multiple  integrations  indicated  in  equation  (5) ,  for 

different  times  tf,  in  order  to  form  the  set  of  overdetermined 

linear  algebraic  equations  whose  solution  gives  the  values 

of  the  parameters  a.,  b.  and  d,  ,  although  only  the  a.  and  b. 

are  of  interest  in  this  work.   The  d,  ,  even  available  are 

k 

not  used. 

The  multiple  integrations  can  be  handled  for  any  numerical 
method  subroutine  and  although  there  are  quite  a  few  highly 
sophisticated  numerical  integration  procedures  available, 
trapezoidal  integration  will  be  used  because  most  of  the 
more  complex  integration  techniques  perform  poorly  when  the 
function  being  integrated  is  discontinuous.   Since  the 
inputs  used  here  are  step-wise  waves,  complex  integration 
techniques  are  not  desirable.   On  the  other  hand,  the  slight 
improvement  in  accuracy  offered  by  more  advanced  methods  is 

not,  in  general,  enough  to  justify  the  tremendous  increase 
in  computational  load  associated  with  their  use. 

The  solution  of  the  set  of  overdetermined  linear  equations 
is  a  classical  problem  in  several  fields  of  mathematics  and 
engineering.   Most  of  the  classical  techniques  for  solving 
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such  problem  are  not  practical.   They  tend  to  magnify 
the  errors  introduced  by  the  finite  precision  of  the  computer 
to  the  point  where  the  solution  is  meaningless.   Fortunately, 
several  modern  methods  are  available  that  display  more 
acceptable  behaviour.   The  method  used  here  was  developed 
by  Golub  [22] .   The  basic  approach  is  to  triangularize  the 
coefficient  matrix  by  performing  a  Choleski  decomposition. 
The  decomposition  is  accomplished  by  applying  repeated 
Householder  transformations  [23]  .   Once  the  coefficient 
matrix  has  been  triangularized  the  unknowns  can  be  obtained 
by  back  substitution.   The  method  is  quite  stable  numerically 
and  is  capable  of  handling  ill-conditioned  coefficient 
matrices. 

B.   CHARACTERISTICS  OF  GOOD . INPUT-OUTPUT  RECORDS 

The  accuracy  with  which  a  system  can  be  identified  is 
strongly  dependent  on  the  input-output  record  used  in  the 
identification.   Since  parameters  are  identified  on  the 
basis  of  their  effect  on  the  output,  it  will  be  impossible 
to  identify  a  parameter  unless  its  effect  is  measurable.   If 
the  input  to  a  system  has  a  frequency  spectrum  that  is  more 
or  less  uniform  over  the  frequency  range  of  interest,  the 
identification  will  probably  be  very  good.   If  the  frequency 
spectrum  of  the  input  is  confined  to  a  narrow  band  of  fre- 
quencies, the  identification  will  probably  be  very  bad.   It 
is  well  known  that  signals  with  sharp  discontinuities  have 
a  broader  bandwidth  than  slowly  varying  continuous  signals. 
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For  this  reason,  the  step-wise  wave  was  chosen  as  the  input 
in  this  research,  after  having  demonstrated  experimentally 
the  preceeding  observations  seeing  that  with  a  step  input 
it  was  never  possible  to  obtain  a  good  identification  even 
for  a  second  order  system.   Periodic  functions  should  also 
be  avoided  at  the  input,  because  very  probably  many  of  the 
algebraic  equations  are  going  to  be  linearly  dependent.   This 
fact  was  also  proved  experimentally;  a  periodic  square  wave 
was  tried  as  input  and  the  identification  was  also  bad  when 
a  step-wise  input  gave  very  good  results. 

Another  desired  characteristics  of  the  input,  experimentally 
proved,  is  that  it  should  be  adjusted  in  such  a  way  that  the 
output,  if  it  is  strictly  increasing,  should  at  some  moment 
in  the  sampling  interval  have  a  negative  second  derivative; 
in  other  words,  the  output  is  desired  to  have  an  inflexion 
point  within  the  period  in  which  the  record  is  taken.   The 
reason  for  this  is  thought  to  lie  in  the  subroutine  DLLSQ 
which  solves  the  overdetermined  set  of  algebraic  equations. 
What  DLLSQ  essentially  does  is  look  for  a  fitting  output 
function  that  minimizes  the  least  square  error  between  the 
generated  and  the  measured  and  if  the  second  and  higher 
order  derivatives  of  the  output  are  positive  at  the  end  of 
the  sampling  interval,  the  function  generated  by  the  sub- 
routine may  in  some  cases  correspond  to  an  unstable  situation. 
Examples  in  Chapter  IV  will  also  demonstrate  this  point. 
In  relation  to  the  sampling  rate  at  which  the  input  and  output 
should  be  digitized,  it  can  be  said  that  low  sampling  rates 
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may  introduce  inaccuracies  in  the  numerical  integrations 
and  as  a  result  poor  identification,  especially  in  the 
location  of  high  frequency  poles  and  zeros.   Maximum 
sampling  rate  available  or  ten  to  one  hundred  times  shorter 
than  the  shortest  system  time  constant  is  recommended.  The 
output  of  the  system  is  also  recommended  to  have  maximum 
value  of  the  order  of  the  maximum  value  of  the  input,  thus 
the  numbers  with  which  the  computer  deals  are  not  very 
different  in  magnitude.   If  the  output  has  not  the  same  order 
of  magnitude  of  the  input,  it  can  as  a  rule  be  amplified 
accordingly,  taking  this  amplification  factor  into  account 
when  the  time  comes  to  see  what  the  actual  identification  is. 

C.   PROGRAMS  DESCRIPTION 

Four  different  computer  programs  were  used  in  this 
research.   A  listing  of  them  is  given  after  Chapter  IV.   The 
first  one  is  the  A/D  program.   This  program  is  used  for  high 
speed  A/D  (analog  to  digital)  conversion.   Its  output  is  a 
seven  track  magnetic  tape  that  contains  the  digitized  input 
and  output  record  of  the  system.   This  program  consists  of 
a  main  driver  program  and  five  subroutines. 

The  source  language  is  FORTRAN  for  the  main  program  and 
METASYMBOL  for  the  subroutines. 

The  A/D  program  provides  a  means  of  sampling  and  digitizing 
analog  information  at  a  rate  (up  to  lOKHz)  determined  by  an 
interrupt  signal  originated  at  the  analog  computer  logic 
patchboard  (Ci  5000) .   (See  Appendix  I) .   The  digital  data 
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are  stored  sequentially  in  one  of  two  buffers  provided  by 
the  main  program.   When  one  buffer  is  full,  the  storage 
sequence  is  continued  at  the  top  of  the  other  buffer.   An 
indicator  notifies  the  main  program  of  this  switching  and 
the  main  program  processes  the  data  in  the  full  buffer 
writing  them  in  the  magnetic  tape  while  the  other  buffer 
is  being  filled. 

The  parameters  to  be  entered  in  the  main  program  are: 

NSAMP  =  Number  of  words  in<  one  block/NCHAN. 

NCHAN  =  Number  of  channels  to  be  digitized. 

NREC   =  Maximum  number  of  records  to  be  digitized  in 
one  file. 

ITAPE  =  Number  of  tape  drive  unit  to  be  used. 

NDEL  =  Delay  between  digital  and  analog  signal  in 
option  7.   The  options  of  this  program  are: 

1  =  Desired  to  input  new  parameters. 

2  =  Start  digitizing  the  analog  input  signals. 

3  =  Write  an  END  OF  FILE  on  the  magnetic  tape. 

4  =  Rewind  the  tape  to  the  loading  point  (beginning) . 

5  =  Skip  the  following  number  of  END  OF  FILE  marks. 

6  =  Print  the  following  amount  of  digitized  data. 

7  =  Actuate  the  A/D  program  for  the  next  single  block 
of  data  encountered  and  feed  the  signal  to  the  recorder. 

The  five  subroutines  are:   ADSTART,  ADSTOP,  ADFAST, 
MTRDY  and  MTOUT. 

Subroutine  ADSTART  initializes  the  A/D  input  operations 
The  calling  sequence  is: 

28 


CALL  ADSTART  (NCH,  BUF ,  NEWBUF ,  NSAM,  RECNUM,  NEXLOC), 
where 

NCH  =  Number  of  A/D  channels  to  sample. 

BUF  =  A  pointer  to  the  first  buffer. 

NEWBUF  =  A  pointer  to  the  second  buffer. 

NSAM  =  Number  of  samples  per  buffer. 

RECNUM  =  A  record  counter  which  is  bumped  by  ADFAST 
each  time  a  buffer  is  filled. 

NEXLOC  =  The  return  address  for  ADFAST  to  be  used  each 
time  a  buffer  is  filled. 

The  first  four  arguments  are  used  to  set  up  pointers 
and  counters  in  ADFAST.   RECNUM  provides  a  means  for  ADFAST 
to  keep  the  main  program  aware  of  the  record  number  of  the 
buffer  currently  being  filled.   NEXLOC  is  an  alternate 
return  location  for  ADFAST.   ADSTART  also  sets  up  the 
interrupt  controller  branch  instructions  and  starts  the  A/D 
hardware.   ADSTOP  stops  the  A/D  operations  and  restores  the 
interrupt  locations  to  their  normal  instructions. 

ADFAST  controls  the  buffer  location  pointer  which  directs 
the  A/D  data  into  core. 

MTRDY  initializes  the  magnetic  tape  handler.   It  has 
the  following  arguments: 

UNIT  =  Magnetic  tape  unit  number. 

BUF1  =  First  word  address  of  first  buffer. 

BUF2  =  First  word  address  of  second  buffer. 

NW   =  Number  of  words  to  be  written  per  record. 

IR   =  Tape  operation  status  word  return. 
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1  if  tape  I/O  channel  is  busy. 

2  if  operation  completed. 

3  if  tape  unit  is  not  ready. 

4  if  tape  writes  error. 

MTOUT  controls  the  stream  of  A/D  data  from  the  buffers 
into  the  tape.   It  is  called  by  the  main  program  each  time 
a  record  of  data  is  ready  for  recording.   The  argument,  NB, 
indicates  which  buffer  is  just  filled  and  is  to  be  output. 
After  setting  the  status  indicator  ilR)    to  1,  to  indicate 
an  operation  not  complete  condition,  the  routine  tests  the 
tape  and  channel  status.   If  either  is  busy,  IR  is  bumped  to 
3  (not  ready  condition)  and  the  routine  exits  back  to  the 
calling  program.   If  everything  is  all  right,  NB  is  checked 
and  the  I/O  instructions  for  writing  the  filled  buffer  on 
tape  are  executed.   If  the  operation  is  completed  without 
error,  IR  is  set  to  2  (operation  completed)  and  to  4  (write 
error)  if  an  error  occurred,  and  execution  is  returned  to 
main  program. 

The  seven  track  magnetic  tape  generated  by  the  A/D 
program  recording  the  digitized  input  and  output  can  not  be 
directly  read  by  the  IBM  360  computer  system,  and  since  this 
system  is  used  to  perform  the  integrations  of  equation  (5) , 
a  seven  to  nine  track  tape  conversion  should  be  done  prior 
to  the  numerical  integration.   This  conversion  is  done  by 
the  IBM  360  using  the  CONVERT  program,  which  is  described 
in  Ref.  21.   This  program  has  several  parameters.   One  of 
them  is  BLKSIZE  which  is  related  with  the  NSAMP  and  NCHAN 
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parameters  used  in  the  A/D  program.   BLKSIZE  can  not  be 
greater  than  32.700  and  for  N€HAN  =  2  limits  the  maximum 
number  of  samples  that  can  be  taken  for  the  input  and  out- 
put record  to  about  4000. 

The  third  program  is  the  main  one  and  is  a  direct 
implementation  of  the  procedure  developed  in  Chapter  II. 
This  program  has  several  important  parameters  that  have  to 
be  defined  at  the  very  beginning.   NP  and  NZ  are  rough 
estimates  of  the  number  of  poles  and  zeros  in  the  system  to 
be  identified.   NP  and  NZ  stand  for  n'  and  m1  in  the  program, 
N'  and  m1  were  defined  in  Chapter  II.   KPTMAX  is  the  number 
of  sample  points  available  in  the  input-output  record.   The 
input  amplitude  is  R  and  the  output  amplitude  is  C.   Because 
the  integration  of  the  time  has  to  be  done,  the  program 
generates  the  time  corresponding  to  each  sample  for  what 
the  sampling  rate  has  to  be  known  and  is  represented  by  the 
parameter  DT  in  the  program.   IPTS  is  the  number  of  sample 
points  that  will  separate  successive  linear  equations. 
Therefore,  the  total  number  of  linear  algebraic  equations 
formed  is  equal  to  KPTMAX/IPTS.   The  program  reads- in  the 
nine  track  tape  generated  by  the  CONVERT  program,  performs 
the  integrations  of  equation  (5)  and  forms  the  overdetermined 
set  of  linear  algebraic  equations.   When  the  set  is  formed, 
subroutine  DLLSQ  is  called.   This  subroutine  is  an  implemen- 
tation of  the  Golub  algorithm  for  solving  overdetermined 
sets  of  linear  equations.   Subroutine  DLLSQ  returns  the 

values  a.,  b.  and  d,  .   The  a.  and  b.  coefficients  are  fed 
l'  3  k        i      j 
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into  the  RTPLSB  subroutine,  which  is  a  polynomial  root 
finder,  in  order  to  obtain  the  poles  and  zeros  of  the  system. 

The  output  of  the  program  are  the  results  of  the 
identification  given  in  both  transfer  function  and  state 
variable  form. 

The  fourth  program  was  written  to  compare  the  outputs 
of  the  actual  system  and  of  the  ones  obtained  from  the 
identification,  when  all  of  them  are  excited  with  a  step. 
This  program  solves  the  differential   equations  that  describe 
the  system.   The  differentials  equations  should  be  arranged 
in  state  variable  form.   The  output  of  the  program  is  a 
drawing   of  the  response  of  the  different  systems  to  be 
compared   in  the  same  plot. 

D.   EXAMPLES 

The  following  examples  show  the  possibilities  of  the 
method,  how  the  order  of  the  system  can  in  general  be 
determined,  if  it  is  not  already  known,  from  a  trial  identi- 
fication run  using  an  estimated  order  greater  than  the 
actual  order  of  the  system.   They  also  show  how  the  accuracy 
of  an  identification  depends  on  the  input  wave  form,  on  the 
R.M.S.  value  of  the  noise  present,  and  as  expected,  does 
not  depend  on  the  initial  conditions.   They  also  indicate 
how  the  method  can  be  used  to  reduce  or  obtain  a  lower  order 
mathematical  model.   In  all  the  examples,  the  MAIN  PROGRAM 
PARAMETERS  were: 

NUMBER  OF  DATA  POINTS     =   4000 
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NUMBER  OF  EQUATIONS       =     40 

DATA  POINTS  PER  EQUATION  =    10  0 

In  examples  1  through  6,  the  input  function  was  a  step- 
wise increasing  wave  form  (staircase  wave  form,  see  Fig.  1, 
Example  1),  with  a  maximum  value  of  about  85  volts.   The 
amplitude  and  the  number  of  steps  were  changed  in  each 
example.   In  Example  3  and  4, the  maximum  amplitude  of  the 
input  was  about  9  volts.   In  examples  6  through  10,  the 
input  function  was  also  a  step-wise  wave  form,  but  instead 
of  being  strictly  increasing,  it  had  one  part  strictly 
increasing  and  the  final  part  strictly  decreasing.  See 
Fig.  4,  Example  5. 

The  systems  investigated  in  Examples  1,  2,  5,  6,  7,  8, 
9,  10  and  11  were  set  up  using  R-C  and  R-C-L  circuits 
isolated  by  operational  amplifiers  which  also  provided  the 
desired  amplification  factor.   In  Examples  3  and  4,  an 
active  system  was  studied.   This  system  was  a  D.C.    positioning 
servo. 

Example  1  was  performed  as  a  first  test  and  in  it  a 
second  order  system  -  two  real  poles  -  was  considered.   The 
input  was  a  step-wise  increasing  function,  (see  Fig.  1), 
with  zero  initial  conditions.   The  results  of  the  identifica- 
tion were  within  the  2%  for  the  dominant  pole  and  the  gain 
constant  and  within  the  10%  for  the  higher  frequency  pole. 
This  example  shows  how  assuming  different  values  for  m'  and 
n' ,  the  results  of  the  identification  gave  a  clear  indication 
of  what  the  actual  order  of  the  system  was. 
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Example  2  illustrates  the  identification  of  a  third 
order  system.   The  input  function  as  in  Example  1  was  a 
step-wise  increasing  function.   The.  initial  conditions  were 
also  zero.   The  identification  of  the  low  frequency  poles 
and  the  gain  constant  were  also  within  3%  of  the  true  values; 
however,  the  pole  at  -100. 0,  was  identified  with  an  error  of 
about  30%,  this  result  was  not  a  surprise  since  the  more 
negative  the  value  of  a  pole,  compared  with  the  other  poles 
in  the  system,  the  less  exact  its  identification.   This 
example  shows  also  the  fact  that  in  some  cases  the  trial 
runs  to  find  out  the  true  values  of  n'  and  m1 ,  did  not  indicate 
very  clearly  what  these  values  were.   In  cases  like  this,  a 
thorough  search  should  be  done  around  the  expected  values  of 
n  and  m  and  different  input-output  records  are  recommended 
to  be  used. 

This  example  also  shows  how  the  method  could  be  used  to 
reduce  the  order  of  the  system,  i.e.,  when  n1  =  3  and  m'  = 
1  were  assumed,  the  results  correspond  practically  to  a 
second  order  system,  since  the  pole  at  -2.22  and  the  zero  at 
-2.20  cancel  each  other.   When  n1  =  2  and  m1  =  0  were 
assumed,  practically  the  same  results  were  obtained.   The 
steady  state  error  coefficient  of  the  second  order  has 
practically  the  same  value  of  the  third  order.   Essentially, 
the  method  neglected  the  effects  of  the  higher  frequency 
pole. 

The  experiment  described  in  Example  3  was  performed 
using  a  D.C.  positioning  servo.   The  transfer  function  of 
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the  servo  was  first  obtained  from  the  asymptotic  approxima- 
tion of  the  closed  loop  frequency  response  plot  shown  in 
Fig.  2.   The  parameters  of  the  servo  were  adjusted  in  such 
a  way  that  two  dominant  complex  roots  were  present.   During 
the  frequency  response  test,  it  was  found  that  the  servo 
had  nonlinearities  since  a  jump  in  the  resonance  frequency 
was  observed  and  also  that  reducing  the  amplitude  of  the 
input  sine  wave  from  twenty  volts  to  one  volt  peak  to  peak 
the  resonance  frequency  was  shifted  from  2.6Hz  to  8.1Hz 
(52  rad/sec) . 

In  taking  the  input-output  record  in  this  example,  the 
input  was  adjusted  to  be  a  step-wise  increasing  function, 
(see  Fig.  1),  but  with  a  step  amplitude  of  about  0.5  volts. 
The  4  volts  step  amplitude  used  in  Examples  1  and  2  was  too 
big;  the  D.C.  motor  ran  continuously. 

The  results  of  the  identification  were  almost  identical 
to  the  ones  obtained  from  the  frequency  response  plot. 

In  Example  4,  the  same  D.C.  servo  was  treated  but  now 
the  parameters  of  the  servo  were  adjusted  in  order  to  have 
two  real  poles.   Also  a  frequency  response  test  was  performed 
and  a  transfer  function  was  obtained  from  the  frequency 
response  curve  asymptotic  approximation  (see  Fig.  3) .   The 
results  of  the  identification  were  very  close  to  the  ones 
found  from  the  frequency  response  plot. 

In  Example  5,  a  third  order  system  is  again  considered, 
but  now  a  zero  was  added. 
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The  input  function  and  the  I.e.  were  the  same  as  in 
Example  1.  The  results  in  this  case  also  showed  a  poor 
identification  of  the  higher  frequency  pole. 

Example  6  shows  the  effect  of  the  input  function  in 
the  accuracy  of  the  identification.   Now  a  fourth  order 
system  is  studied  using  two  different  inputs.   The  first 
input  was  the  same  as  in  Example  1  and  is  shown  together 
with  the  output  of  the  system  in  Fig.  4.   With  this  input, 
the  results  of  the  identification  were  very  poor  since  not 
a  single  pole  was  correctly  identified.   Several  fourth 
order  systems  were  treated  with  this  type  of  input  and  the 
identification  was  bad,  giving   in  some  cases  poles  in  the 
right  half  plane.   It  was  observed,  however,  that  the  output 
of  the  system,  at  the  end  of  the  record  had  a  positive 
second  derivative  (see  Fig.  4) .   It  was  thought  that  perhaps 
this  fact  could  be  the  cause  of  the  poor  identification 
considering  that  an  unstable  system  might  have  the  same 
output,  for  the  finite  interval  of  time  considered,  and  this 
unstable  system  could  be  the  one  identified  or  generated 
by  the  DLLSQ  subroutine.   With  the  preceding  in  mind,  the 
input  function  was  changed  and  adjusted  in  such  a  way  that 
the  output  had  the  second  derivative  negative  at  a  certain 
instant  within  the  time  interval  considered  (see  Fig.  5). 
The  results  in  this  case  were  very  good  since  even  the  high 
frequency  pole  was  identified  with  an  error  within  the  6% 
of  the  true  value.   The  I.e.  were  also  zero  value. 


36 


In  Example  7 ,  the  identification  of  a  system  with  four 
poles  and  two  zeroes  is  presented.  The  input  function  was 
as  the  one  shown  in  Fig.  5.  The  results  in  this  case  were 
also  very  good. 

Example  8  shows  how  the  I.C.  does  not  affect  the  identi- 
fication, the  results  being  practically  the  same  with  I.C. 
equal  to  zero  as  with  I.C.  0;  several  different  initial 
conditions  were  considered. 

In  Example  9,  a  third  order  system  was  studied,  but  noise 
was  added  at  the  input  and  output.   The  results  showed  the 
great  influence  of  the  noise  in  the  identification.   It  was 
found  that  in  order  to  have  a  good  identification,  the 
maximum  R.M.S.  value  of  the  nois^f  should  not  be  greater  than 
about  .5  volts.   1.0  volts  R.M.S.  noise  made  the  identifica- 
tion very  bad.   The  noise  added  was  a  gaussian,  low  frequency 
noise  -  from  a  few  cycles  per  second  to  about  300  cycles  per 
second  -  with„z,ero  mean. 

Examples  10  and  11  show  how  the  method  can  be  used  to 
reduce  the  order  of  the  system.   In  Example  10,  a  fifth 
order  system  was  treated.   The  identification  was  not  good 
for  the  higher  frequency  poles.   Other  fifth  order  systems 
were  studied  and  in  some  cases  none  of  the  poles  was  identi- 
fied, which  indicates  that  fourth  order  systems  are  perhaps 
the  highest  order  systems  that  can  be  identified  with  the 
kind  of  inputs  used  in  this  work. 

In  Example  10,  the  reduction  was  not, as  a  direct  conse- 
quence of  the  poor  identification,  very  good  as  can  be  seen 
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comparing  the  responses  to  a  step  input  of  the  different 
models,  see  Fig.  6.   However,  it  is  observed  that  the 
three  models  generated  by  the  method,  5,3,  and  2 
order,  give  responses  fairly  close  to  one  another. 

In  Example  11,  although  the  identification  is  considered 
very  good,  the  responses  of  the  different  reduced  models 
practically  coincide  during  the  transient  period,  but  in 
steady  state  an  error  of  about  30%  is  observed  between  the 
step  responses  of  the  fourth  order  identified  model  and  the 
third  and  second  order  reduced  models  (see  Fig.  7).   However, 
if  the  gain  constant  of  the  reduced  models  is  adjusted  in 
order  to  have  the  same  steady  state  gain  as  the  fourth  order 
model,  the  responses  are  practically  the  same  (see  Fig.  8). 

It  is  to  be  noticed  that  the  systems  considered  true, 
in  the  examples,  were  true  within  the  tolerances  of  the  values 
of  the  resistors,  capacitors,  coils  and  amplifiers  used. 
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EXAMPLE  I 


SYSTEM  TO  BE  IDENTIFIED 
SYSTEM  TRANSFER  FUNCTION 

POLES  REAL 

1  -10.0 

2  -100.0 


ZEROES 


REAL 


IMAGINARY 

0.0        J 
0.0        J 

IMAGINARY 


GAIN  CONSTANT  =  1000.0 


TESTS  FOR  EXAMPLE  I 


RUN 

NO. 

OF 

POLES 

NO. 

OF  ZEROES 

RESULTS 

a 

2 

0 

Good 

b 

3 

■ 

1 

See  Table  of  Data 

c 

4 

2 

See  Table  of  Data 

d 

2 

1 

See  Table  of  Data 

INPUT  FUNCTION  STEP-WISE  STRICTLY  INCREASING  SHOWN  IN 

FIG.  1. 

INITIAL  CONDITIONS  EQUAL  ZERO. 
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INPUT 
VOLTS 

80     , 


4  00    MS 


rt 


I2TPUT  FU1TCTION  USZD   III  EXAMPLE  ::0;  1 


OUTPUT 
VOLTS 


OUTPUT  FUNCTION   EXAI-IPLE  ITO.l. 


FIGURE   1 


4  00  MS 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -9.8488868 

2  -111.04012 


IMAGINARY 
0.0 
0.0 


J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  984.74146 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 

1  1093.6215 

2  120.88900 


B  VECTOR 
2 


984.74146 


C  VECTOR 

1 


1.0000000 


REMARKS : 

ASK  FOR  TWO  POLES  NO  ZEROES. 
GOOD  RESULTS. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (b) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-97.994740 

0.0 

J 

2 

2.2362954 

0.0 

J 

3 

-9.8646060 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

2.2427425 

0.0 

J 

GAIN  CONSTANT  =  867.90796 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 

1  -2161.7809 

2  725.47414 

3  105.62305 


B  VECTOR 

3 

867.90796 

C  VECTOR 

1 

-2.2427425 

2 

1.0000000 

REMARKS : 

ASK  FOR  THREE  POLES  ONE  ZERO. 

RESULTS  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE  SYSTEM 

SINCE  ZERO  AT  2.24  CANCELS  POLE  AT  2.23. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (c) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-4.4110354 

-54.484437 

J 

2 

-4.4110354 

54.484437 

J 

3 

-99.827986 

0.0 

J 

4 

-9.8612480 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

-7.0865674 

-59.922949 

J 

2 

-7.0865674 

59.922949 

J 

GAIN  CONSTANT  =  727.27368 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

2941483.4 

2 

336437.35 

3 

4940.1259 

4 

118.51130 

B 

VECTOR 

4 

727.27368 

C 

VECTOR 

1 

3640.9793 

2 

14.173135 

3 

1.0000000 
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REMARKS : 

ASK  FOR  4  POLES  TWO  ZEROES. 

RESULTS  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE 
SYSTEM,  SINCE  COMPLEZ  ZEROES  CANCEL  COMPLEX  POLES. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (d) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-9.8592476 

0.0 

J 

2 

-97.619069 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

-847.73348 

0.0 

J 

GAIN  CONSTANT  =  1.0221453 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

962.47832 

2 

107.47832  ' 

B 

VECTOR 

2 

1.0221463 

C 

VECTOR 

1 

847.73348 

2 

1.0000000 

REMARKS : 

ASK  FOR  TWO  POLES  ONE  ZERO. 

RESULTS  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE  SYSTEM. 
SINCE  THE  ZERO  AT  -  847.7  CAN  BE  NEGLECTED. 
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EXAMPLE  2 

SYSTEM  TO  BE  IDENTIFIED 

SYSTEM  TRANSFER  FUNCTION 

POLES             REAL  IMAGINARY 

1  -1.000  0.0          J 

2  -10.00  0.0          J 

3  -100.0  0.0          J 

ZEROES 

GAIN  CONSTANT  -  10000.0 

TESTS  FOR  EXAMPLE  2 


RUN 

NO. 

OF 

POLES 

NO. 

OF 

ZEROES 

RESULTS 

a 

3 

0 

Good 

b 

4 

1 

See 

Table 

of 

Data 

c 

4 

0 

See 

Table 

of 

Data 

d 

4 

2 

See 

Table 

of 

Data 

e 

3 

1 

See 

Table 

of 

Data 

f 

3 

1 

See 

Table 

of 

Data 

INPUT  FUNCTION  STEP-WISE  STRICTLY  INCREASED  SIMILAR 
TO  THE  ONE  SHOWN  IN  FIG.  1. 

INITIAL  CONDITIONS  EQUAL  ZERO. 
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IDENTIFICATION  OF  UNKNOV7N  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -1.1072321 

2  -10.224865 

3  -131.11147 


IMAGINARY 
0.0 
0.0 
0.0 


J 
J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  11153.488 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

1484.3522 

2 

1497.0893 

3 

142.44357 

B  VECTOR 

3 

11153.488 

C  VECTOR 

1.0000000 

REMARKS : 

ASK  FOR  THREE  POLES  NO  ZEROES. 

GOOD  RESULTS  FOR  LOW  FREQUENCY  POLES,  NOT  SO  GOOD  FOR 
HIGH  FREQUENCY  POLE. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (b) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-1.1675351 

0.0 

J 

2 

-9.9863699 

0.0 

J 

3 

-1.7691065 

-59.923677 

J 

ZEROES 

REAL 

IMAGINARY 

1 

157.78387 

0.0 

J 

GAIN  CONSTANT  =  -1919.7240 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

41903.747 

2 

40128.130 

3 

3645.1011 

4 

14.692118 

B 

VECTOR 

4 

-1919.1240 

C 

VECTOR 

1 

-157.78387 

2 

1.0000000 

REMARKS : 

ASK  FOR  FOUR  POLES  ONE  ZERO. 

RESULTS  TO  NOT  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE 
SYSTEM,  INDICATING  FOUR  POLES  WHEN  THERE  ARE  THREE.   THE 
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ZERO  IN  THE  RIGHT  HALF  PLANE  SHOULD  BE  NEGLECTED 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (c) 


SYSTEM  TRANSFER  FUNCTION 
POLES        .     REAL 

1  -1.1300231 

2  -10.158884 

3  -11.337775 

4  -11.337775 


IMAGINARY 

0.0 

0.0 
-53.285012 
+53.285012 


J 
J 
J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  252014.00 

SYSTEM  STATE  VARIABLE  (PHASE  FORM) 
A  VECTOR 


1 

34070.104 

2 

33763.954 

3 

3235.2996 

4 

33.964457 

B  VECTOR 

4 

252014.00 

C  VECTOR 

1 

1.0000000 

REMARKS : 

ASK  FOR  FOUR  POLES  NO  ZEROES. 

RESULTS  DO  NOT  GIVE  CLEAR  INDICATION  OF  ORDER  OF  THE 
SYSTEM  INDICATING  FOUR  POLES  WHEN  THERE  ARE  THREE. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (d) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-1.1196795 

0.0 

2 

-10.190750 

0.0 

3 

0.16021282 

-59.917890 

4 

0.16021282 

59.917890 

ZEROES 

REAL 

IMAGINARY 

1 

3.5557973 

-30.108372 

2 

3.5557973 

30.108372 

J 
J 
J 
J 


J 
J 


GAIN  CONSTANT  =  332.37549 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

40965.288 

2 

40602.813 

3 

3597.9654 

4 

10.990004 

B 

VECTOR 

4 

332.37549 

C 

VECTOR 

1 

919.15775 

2 

-7.1115946 

3 

1.00000000 
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REMARKS : 

ASK  FOR  FOUR  POLES  TWO  ZEROES 

RESULTS  INDICATE  TWO  POLES  WHEN  THERE  ARE  THREE 
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IDENTIFICATION  OF  UNKNOV7N  SYSTEM  -  Run  (e) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-1.0459484 

0.0 

J 

2 

-10.251872 

0.0 

J 

3 

2.2244020 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

2.2033545 

0.0 

J 

GAIN  CONSTANT  =  83.810272 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

23.852104 

2 

35.853822 

3 

13.522222 

B  VECTOR 

3 

83.810272 

C  VECTOR 

1 

2.2033545 

2 

1.0000000 

REMARKS : 

ASK  FOR  THREE  POLES  ONE  ZERO. 

RESULTS  INDICATE  TWO  POLES  WHEN  THERE  ARE  THREE. 
POLE  AT  -2.2  IS  CANCELLED  BY  ZERO  AT  -2.22. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (f) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -1.0934113 

2  -10.207312 


IMAGINARY 
0.0 
0.0 


J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  90.3204350 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 

1  11.232745 

2  11.300723 
B  VECTOR 

2  90.320435 

C  VECTOR 

1  1.0000000 

REMARKS : 

ASK  FOR  TWO  POLES  NO  ZEROES. 

RESULTS  INDICATE  TWO  POLES  NO  ZEROES  WHEN  THERE  ARE 
THREE.   METHOD  NEGLECTED  POLE  AT  -100.0,  OTHERWISE  IDENTI- 
FICATION VERY  GOOD. 
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EXAMPLE  3 


SYSTEM  TO  BE  IDENTIFIED 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -7.8 

2  -7.8 
ZEROES  REAL 


IMAGINARY 
...-51.  5 
51.5 
IMAGINARY 


J 
J 


GAIN  CONSTANT  =  295  0.0 


TESTS  FOR  EXAMPLE  3 


RUN 

NO. 

OF  POLES 

NO. 

OF  RESULTS 

RESULTS 

a 

2 

0 

Good 

D.C.  POSITIONING  SERVO  FREQUENCY  RESPONSE  IN  FIG.  2 
INPUT  FUNCTION  STEP-WISE  INCREASING  SIMILAR  TO  THE  ONE 
SHOWN  IN  FIG.  1,  BUT  WITH  STEP  AMPLITUDE  ABOUT  0.5  VOLTS 

INITIAL  CONDITIONS  EQUAL  ZERO. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -5.3585933 

2  -5.3585933 


IMAGINARY 

-56.771298 

56.771298 


J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  3332.2095 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 

1  3251.6948 

2  10.717187 
B  VECTOR 

2  3332.2095 

C  VECTOR 

1  1.0000000 

REMARKS: 

ASK  FOR  TWO  POLES  NO  ZEROES. 
GOOD  RESULTS. 


57 


EXAMPLE  4 


SYSTEM  TO  BE  IDENTIFIED 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -4.9 

2  -32.0 


IMAGINARY 
0.0 
0.0 


J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  150.0 


TESTS  FOR  EXAMPLE  4 


RUN      NO.  OF  POLES 
a  2 


NO.  OF  ZEROES 
0 


RESULTS 
Good 


D.C.  POSITIONED  SERVO. 

FREQUENCY  RESPONSE  IN  FIG.  3. 

INPUT  FUNCTION  THE  SAME  AS  IN  EXAMPLE  2. 

INITIAL  CONDITIONS  EQUAL  ZERO. 
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IDENTIFICATION  OF  UNKNOV7N  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -5.2129680 

2  -36.525846 


IMAGINARY 
0.0 
0.0 


J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  178.37120 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 

1  190.40807 

2  41.738815 
B  VECTOR 

2  178.37120 

C  VECTOR 

1  1.0000000 

REMARKS : 

ASK  FOR  TWO  POLES  NO  ZEROES. 
GOOD  RESULTS. 
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EXAMPLE  5 


SYSTEM  TO  BE  IDENTIFIED 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -1.0 

2  -10.0 

3  -55.0 


ZEROES 

1 


REAL 
-5.0 


GINARY 

0.0 

J 

0.0 

J 

0.0 

J 

GINARY 
0.0 

J 

GAIN  CONSTANT  =  500.0 


TESTS  FOR  EXAMPLE  5 


RUN      NO.  OF  POLES 
a  3 

b  2 


NO.  OF  ZEROES      RESULTS 
1  Good 

0       See  Table  of  Data 


INPUT  FUNCTION  STEP-WISE  STRICTLY  INCREASING  SIMILAR  TO 
THE  ONE  SHOWN  IN  FIG.  1. 
INITIAL  CONDITIONS  EQUAL  ZERO. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-0.89534766 

0.0 

J 

2 

-8.3914174 

0.0 

J 

3 

-88.600646 

p.o 

J 

ZEROES 

REAL 

IMAGINARY 

1 

-3.8670622 

0.0 

J 

GAIN  CONSTANT  =  706.88770 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

665.67756' 

2 

830.32662 

3 

97.887411 

B  VECTOR 

3 

706.88770 

C  VECTOR 

1 

3.8670622 

2 

1.0000000 

REMARKS : 

ASK  FOR  THREE  POLES  ONE  ZERO. 

RESULTS  SHOW  POOR  IDENTIFICATION  OF  THE  HIGHER 

FREQUENCY  POLE. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (b) 


SYSTEM  TRANSFER  FUNCTION 

POLES  REAL  IMAGINARY 

1  -3.2441577         0.0 

2  -108.41738         0.0 


J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  763.43213 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 

1  351.72309 

2  111.66154 
B  VECTOR 

2  763.43213 

C  VECTOR 

1  1.0000000 

REMARKS : 

ASK  FOR  TWO  POLES  NO  ZEROES. 

RESULTS  CAN  BE  CONSIDERED  AS  A  LOWER  ORDER  MATHEMATICAL 
MODEL. 


63 


EXAMPLE  6 


SYSTEM  TO  BE  IDENTIFIED 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -2.0 

2  -10.0 

3  -33.0 

4  -100.0 


ZEROES 


REAL 


IMAGINARY 

0.0 

J 

0.0 

J 

0.0 

J 

0.0 

J 

IMAGINARY 

GAIN  CONSTANT  =  66  0000 


TESTS  FOR  EXAMPLE  6 


RUN      NO.  OF  POLES 

a  4 

b  4 

c  6 

d  4 

e  5 

f  4 

g  5 

h  3 


NO.  OF  ZEROES 

0 
0 
2 
2 
2 
1 
1 
0 


RESULTS 

Bad 
Good 
See  Table  of  Data 
See  Table  of  Data 
See  Table  of  Data 
See  Table  of  Data 
See  Table  of  Data 
See  Table  of  Data 


INPUT  FUNCTION  FOR  RUN  (a)  STEP-WISE  STRICTLY  INCREASING 
SHOWN  IN  FIG.  4. 

INPUT  FUNCTION  FOR  RUNS  (b) ,  (c) ,  (d) ,  (e) ,  (f )  ,  (g)  ,  and 
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(h)  STEP-WISE  INCREASING-DECREASING  SHOWN  IN  FIG.  5 
INITIAL  CONCLUSION  EQUAL  ZERO. 
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INPUT 
VOLTS 

80 


S 


S 


TIME  -  MSEC 

FIRST  INPUT  FUNCTION  USED  IN  EXAMPLE  NO.   6 


400   • 


OUTPUT 
VOLTS 


TIME    -    MSEC 


OUTPUT   COHTt^POHLSNT  TO    IKA  ABOViU   INPUT- 


FIGURE 


400 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -1.5013900 

2  -1.5013900 

3  -5.4821332 

4  -5.4321332 


ZEROES 


REAL 


IMAGINARY 

-4.3496191 

J 

4.3496191 

J 

-26.089080 

J 

26.089080 

J 

IMAGINARY 

GAIN  CONSTANT  =  62546.4  30 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

•  15047.775 

2 

2366.2077 

3 

764.79050 

4 

13.967046 

B  VECTOR 

4 

62546.430 

C  VECTOR 

1 

1.0000000 

REMARKS : 

INPUT  FUNCTION  -  STEP-WISE  INCREASING  FUNCTION  SHOWS  IN 
FIG.  4. 

ASK  FOR  FOUR  POLES  NO  ZEROES. 
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RESULTS  NOT  GOOD. 

NOTE  THE  POSITIVE  SECOND  DERIVATIVE  OF  THE  OUTPUT  AT 
THE  END  OF  THE  SAMPLING  INTERVAL. 
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INPUT 
VOLTS 

80 


TIME  -  MSEC 


400 


OUTPUT 
VOLTS 


80 


SECOND  INPUT  FUITCTIGIT  USED   III  EXAMPLE  ITO.      6 


TIME    -    MSEC 


400 


OUTPUT   CORRESPOITDSITT  TO 


above  i:;put 


FIGURE 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (b) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -2.0172066 

2  -10.764637 

3  -33.596973 
..  4  -105.65671 


ZEROES 


REAL 


IMAGINARY 

0.0 

J 

0.0 

J 

0.0 

J 

0.0 

J 

IMAGINARY 

GAIN  CONSTANT  =  612357.94 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

77080.941 

2 

48396.117 

3 

5351.3789 

4 

152.03553 

B  VECTOR 

4 

612357.94 

C  VECTOR 

1 

1.0000000 

REMARKS : 

INPUT  FUNCTION  STEP-WISE  INCREASING-DECREASING  SHOWN 
IN  FIG.  5. 

ASK  FOR  FOUR  POLES  NO  ZEROES. 
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RESULTS  VERY  GOOD. 

NOTE  THE  NEGATIVE  SECOND  DERIVATIVE  DURING  SOME  TIME 
WITHIN  THE  SAMPLING  INTERVAL. 
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IDENTIFICATION  OF  UNKNOV7N  SYSTEM  -  Run  (c) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-2.0375228 

0.0 

J 

2 

-10.649015 

0.0 

J 

3 

-34.253558 

0.0 

J 

4 

-101.99791 

0.0 

J 

5 

6.7412854 

-96.372163 

J 

6 

6.7412854 

96.372163 

J 

ZEROES 

REAL 

IMAGINARY 

1 

6.2635645 

-96.069197 

J 

2 

6.2635645 

-96.069197 

J 

GAIN  CONSTANT  =  603327.00 


SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 

0.7075898D   09 
0.44024818D   09 
48381252.0 


1 
2 
3 
4 
5 
6 
B  VECTOR 
6 


1366621.4 
12569.020 

135.45544 

603327.00 
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* 


C  VECTOR 

1 

9268.5229 

2 

-12.527129 

3 

1.0000000 

REMARKS: 

ASK  SIX  POLES  TWO  ZEROES. 

RESULTS  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE 
SYSTEM,  SINCE  COMPLEX  POLES  ARE  CANCELLED  BY  COMPLEX  ERRORS 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (d) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-2.0305074 

0.0 

J 

2 

-10.677930 

0.0 

J 

3 

-34.308456 

0.0 

J 

4 

-96.072665 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

-516.60587 

0.0 

J 

2 

751.83224 

0.0 

J 

GAIN  CONSTANT  =  -1.4570827 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

71464.878 

2 

44715.215 

3 

4974.7267 

4 

143.08956 

B  VECTOR 

4 

-1.4570827 

C  VECTOR 

1 

-388400.95 

2 

-235.22636 

3 

1.0000000 

74 


REMARKS : 

ASK  FOR  FOUR  POLES  TWO  ZEROES. 

RESULTS  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE 
SYSTEM,  SINCE  THE  TWO  EXCESS  ZEROES  APPEAR  ONE  IN  THE 
R.H.P.  AND  THE  OTHER  IN  THE  L.H.P.  BUT  CAN  BE  NEGLECTED. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (e) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-2.7344050 

0.0 

J 

2 

-8.8026617 

0.0 

J 

3 

-43.843968 

0.0 

J 

4 

-68.587703 

0.0 

J 

5 

8.0082667 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

7.9267041 

0.0 

J 

2 

-622.31546 

0.0 

J 

GAIN  CONSTANT  =  777.31934 


SYSTEM  STATE 

VARIABLES  (PHASE  FORM) 

A  VECTOR 

1 

-579657.55 

2 

-227126.83 

3 

2737.3550 

4 

3335.5841 

5 

115.96047 

B  VECTOR 

5 

777.31934 
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REMARKS ; 

ASK  FOR  FIVE  POLES  TWO  ZEROES 

RESULTS  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE 
SYSTEM,  SINCE  THE  EXCESS  POLE  AT  8.0  IS  CANCELLED  BY  THE 
ZERO  AT  +7.9  AND  THE  OTHER  EXCESS  ZERO  AT  -  622.3  CAN  BE 
NEGLECTED. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (f) 


SYSTEM  TRANSFER  FUNCTIONS 


POLES 

REAL 

IMAGINARY 

1 

-2.0267057 

0.0 

J 

2 

-10.704298 

0.0 

J 

3 

-34.037221 

0.0 

J 

4 

-101.07630 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

-4045.2488 

0.0 

J 

GAIN  CONSTANT  =  146.26987 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

74636.682 

2 

46730.407 

3 

5182.1817 

4 

147.84453 

B  VECTOR 

4 

146.26987 

C  VECTOR 

1 

4045.2488 

2 

1.0000000 

REMARKS : 

ASK  FOR  FOUR  POLES  ONE  ZERO 

RESULTS  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE 
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SYSTEM,  SINCE  THE  EXCESS  ZERO  ZT  -  4045.2  CAN  BE  NEGLECTED 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (g) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-2.1810404 

0.0 

J 

2 

-10.145443 

0.0 

J 

3 

-35.910038 

0.0 

J 

4 

-96.915087 

0.0 

J 

5 

16.571798 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

16.789359 

0.0 

J 

GAIN  CONSTANT  =  580130.69 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

-1276179.1 

2 

-682609.54 

3 

-39334.691 

4 

2734.1957 

5 

128.57981 

B  VECTOR 

5 

580130.69 

C  VECTOR 

1 

-16.789359 

2 

1.0000000 
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REMARKS : 

ASK  FOR  FIVE  POLES  ONE  ZERO. 

RESULTS  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE 
SYSTEM,  SINCE  EXCESS  POLE  AT  16.5  IS  CANCELLED  BY  EXCESS 
ZERO  AT  16.7. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (h) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -1.0271020 

2  -15.051165 

3  -15.051165 


IMAGINARY 
0.0 
-9.0280041 
9.0280041 


J 

J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  4051.1528 


SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 
2 
3 

B  VECTOR 

3 
C  VECTOR 

1 


316.39101 
338.96060 
31.129433 

4051.1528 

1.0000000 


REMARKS : 

ASK  FOR  THREE  POLES  NO  ZEROES. 

RESULTS  CAN  BE  CONSIDERED  AS  A  LOWER  ORDER  MATHEMATICAL 
MODEL. 
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EXAMPLE  7 


SYSTEM  TO  BE  IDENTIFIED 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -25.0 

2  -67.0 

3  -150 

4  -150 


ZEROES 
1 
2 


REAL 

-3.1 

-10.0 


IMAGINARY 

0.0 

J 

0.0 

J 

-185.3 

J 

185.3 

J 

IMAGINARY 

0.0 

J 

0.0 

J 

TESTS  FOR  EXAMPLE  7 


RUN 

NO. 

OF  POLES 

a 

6 

b 

4 

c 

4 

NO.  OF  ZEROES     RESULTS 

3       See  Table  of  Data 

2  Good 

3  See  Table  of  Data 


INPUT  FUNCTION  SIMILAR  TO  THE  ONE  SHOWN  IN  FIG.  5. 
INITIAL  CONDITIONS  EQUAL  TO  ZERO. 
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IDENTIFICATION  OF  UNKNOV7N  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-26. 773050 

0.0 

J 

2 

-72.749634 

0.0 

J 

3 

-101.14084 

173.81324 

J 

4 

-101.14084 

-173.81324 

J 

5 

960.71376 

0.0 

J 

6 

52.080470 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

50.693422 

0.0 

J 

2 

-3.0909162 

0.0 

J 

3 

-9.5071748' 

0.0 

J 

GAIN  CONSTANT  =  -0.33196749E   09 


SYSTEM  STATE 

VARIABLES  (PHASE 

FORM) 

A  VECTOR 

1 

0.39410705D 

13 

2 

0.14131408D 

12 

3 

-0.12683604D 

10 

4 

-43800405. 

5 

-193111.44 

6 

-710.98986 

B  VECTOR 

6 

-0.33196749E 

09 

84 


C  VECTOR 

1 

-1489.6708 

2 

-609.25446 

3 

-38.095331 

4 

1.0000000 

REMARKS : 

ASK  FOR  SIX  POLES  THREE  ZEROES. 

RESULTS  INDICATE  FOUR  POLES  TWO  ZEROES,  SINCE  POLE 
AT  52.0  IS  CANCELLED  BY  ZERO  AT  50.7  AND  POLE  AT  960.7 
CAN  BE  NEGLECTED. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (b) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-26.992965 

0.0 

J 

2 

-77.607958 

0.0 

J 

3 

-126.52085 

-182.98579 

J 

4 

-126.52085 

182.98579 

J 

ZEROES 

REAL 

IMAGINARY 

1 

-3.1577742 

2 

-9.376548. 

0.0 

J 

GAIN  CONSTANT  =  444891.63 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

0.10367784D 

2 

5706927.6 

3 

78054.592 

4 

357.64263 

B  VECTOR 

4 

444891.63 

C  VECTOR 

1 

29.609022 

2 

12.534322 

3 

1.0000000 
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REMARKS : 

ASK  FOR  FOUR  POLES  TWO  ZEROES. 
GOOD  RESULTS. 
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IDENTIFICATION  OF  UNKNOV7N  SYSTEM  -  Run  (c) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

•  1 

-27.385227 

0.0 

J 

2 

-75.010363 

0.0 

J 

3 

-119.37880 

-181.29378 

J 

4 

-119.37880 

181.29378 

J 

ZEROES 

REAL 

IMAGINARY 

1 

-3.1367869 

0.0 

J 

2 

-9.4693811 

0.0 

J 

3 

-5036.4256 

0.0 

J 

GAIN  CONSTANT  =  82.045258 


SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 
2 
3 

4 
B  VECTOR 

4 
C  VECTOR 

1 

2 

3 


96790162.0 
5315200.5 
73620.633 
341.15318 

82.045258 

149599.12 
63519.731 
5049.0317 
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4  1.0000000 

REMARKS : 

ASK  FOR  FOUR  POLES  THREE  ZEROES. 

RESULTS  GIVE  CLEAR  INDICATION  OF  THE  ORDER  OF  THE 
SYSTEM,  SINCE  ZERO  AT  -5036.4  CAN  BE  NEGLECTED. 


89 


EXAMPLE  8 


SYSTEM  TO  BE  IDENTIFIED 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -1.0 

2  -10.0 

3  25.0 


IMAGINARY 
0.0 
0.0 
0.0 


J 
J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  2  50  0 


TESTS  FOR  EXAMPLE  8 


RUN 

NO. 

OF  POLES 

NO. 

OF  ZEROES 

RESULTS 

a 

3 

0 

Good 

b 

3 

0 

Good 

c 

3 

0 

Good 

INPUT  FUNCTION  SIMILAR  TO  THE  ONE  SHOWN  IN  FIG.  5. 
INITIAL  CONDITIONS  DIFFERENT  FROM  ZERO  FOR  THE  THREE 
DIFFERENT  RUNS. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

•  1  -0.97837856 

2  -10.976897 

3  -24.694433 


ZEROES 


REAL 


IMAGINARY 

0.0 

J 

0.0 

J 

0.0 

J 

IMAGINARY 

GAIN  CONSTANT  =  2  378.9880 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

2 

3 
B  VECTOR 

3 
C  VECTOR 

1 


265.20737 
305.93833 
36.649709 

2378.9880 

1.0000000 


REMARKS: 

ASK  FOR  THREE  POLES  NO  ZEROES. 

GOOD  RESULTS  WITH  I.C.  DIFFERENT  FROM  ZERO 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (b) 


SYSTEM  TRANSFER  FUNCTION 

POLES  REAL  IMAGINARY 

1  -0.96861216         0.0 

2  -10.956655         0.0 

3  -24.782535         0.0 


J 
J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  2379.6726 


SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 
2 
3 

B  VECTOR 

3 
C  VECTOR 

1 


263.01082 

306.15108' 

36.707801 

2379.6726 

1.0000000 


REMARKS : 

ASK  FOR  THREE  POLES  NO  ZEROES. 

GOOD  RESULTS  WITH  I.C.  DIFFERENT  FROM  ZERO 
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IDENTIFICATION  OF  UNKNOV7N  SYSTEM  -  Run  (c) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -24.839515 

2  -0.95989923 

3  -10.973508 


IMAGINARY 
0.0 
0.0 
0.0 


J 
J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  2386.2974 


SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 
2 
3 

B  VECTOR 

3 
C  VECTOR 

1 


261.64609 
306.95352 
36.772923 

2386.2974 

1.0000000 


REMARKS: 

ASK  FOR  THREE  POLES  NO  ZEROES. 

GOOD  RESULTS  WITH  I.C.  DIFFERENT  FROM  ZERO 
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EXAMPLE  9 


SYSTEM  TO  BE  IDENTIFIED 


SAME  SYSTEM  EXAMPLE  8. 


TESTS  FOR  EXAMPLE  9 


RUN 

NO. 

OF  POLES 

NO. 

OF  ZEROES 

RESULTS 

a 

3 

0 

Bad 

b 

3 

0 

Almost  Good 

c 

3 

0 

Good 

d 

3 

0 

Good 

INPUT  FUNCTION  STEP-WISE  INCREASING-DECREASING  SIMILAR 
TO  THE  ONE  SHOWN  IN  FIG.  5. 

INITIAL  CONDITIONS  EQUAL  ZERO. 

LOW  FREQUENCY  NOISE  -  FROM  A  FEW  CYCLES  TO  ABOUT  300. 
CYCLES  PER  SECOND  -  WAS  ADDED  AT  THE  INPUT  AND  OUTPUT. 

DIFFERENT  R.M.S.  VALUES  OF  NOISE  WERE  CONSIDERED. 
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IDENTIFICATION  OF  UNKNOV7N  SYSTEM  -  Run  (a) 

SYSTEM  TRANSFER  FUNCTION 

POLES  REAL  IMAGINARY 

1  -0.25785674D-01       0.0  J 

2  -16.598045   -7.2418627       J 

3  -16.598045    7.2418627       J 

ZEROES  REAL  IMAGINARY 

GAIN  CONSTANT  =  2294.0054 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

8.4561459 

2 

328.79567 

3 

33.221877 

B  VECTOR 

3 

2294.0054 

C  VECTOR 

1 

1.0000000 

REMARKS : 

ASK  FOR  THREE  POLES  NO  ZEROES. 
R.M.S.  VALUE  OF  THE  NOISE  1.0  VOLTS. 
RESULTS  NOT  GOOD. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  RUN  (b) 


SYSTEM  TRANSFER  FUNCTION 
.   POLES  REAL 

1  -23.206882 

2  -12.305073 

3  ,        -0.72464570 


ZEROES 


REAL 


IMAGINARY 

0.0 

J 

0.0 

J 

0.0 

J 

IMAGINARY 

GAIN  CONSTANT  =  2364.4692 


SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 
2 

3 

B  VECTOR 

3 
C  VECTOR 

1 


206.93154 

311.29595' 

36.236600 

2364.4692 

1.0000000 


REMARKS : 

ASK  FOR  THREE  POLES  NO  ZEROES. 
R.M.S.  VALUE  OF  THE  NOISE  .75  VOLTS 
RESULTS  ALMOST  GOOD. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (c) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -1.1475703 

2  -10.808421 

3  -23.475558 


ZEROES 


REAL 


IMAGINARY 

0.0 

J 

0.0 

J 

0.0 

J 

IMAGINARY 

GAIN  CONSTANT  =  2289.9822 


SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 
2 
3 

B  VECTOR 

3 
C  VECTOR 

1 


291.17727 
293.07699 
35.431549 

2289.9822 

1.0000000 


REMARKS : 

ASK  FOR  THREE  POLES  NO  ZEROES. 
R.M.S.  VALUE  OF  THE  NOISE  0.5  VOLTS 
RESULTS  GOOD. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (d) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -1.0495050 

2  -10.680594 

3  -25.023091 


ZEROES 


REAL 


IMAGINARY 

0.0 

J 

0.0 

J 

0.0 

J 

IMAGINARY 

GAIN  CONSTANT  =  2381.8215 


SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 
2 
3 

B  VECTOR 

3 
C  VECTOR 

1 


280.49226 
304.73267 
36.753190 

2381.8215 

1.0000000 


REMARKS : 

ASK  FOR  THREE  POLES  NO  ZEROES. 
R.M.S.  VALUE  OF  THE  NOISE  .25  VOLTS 
RESULTS  VERY  GOOD. 


98 


EXAMPLE  10 


SYSTEM  TO  BE  IDENTIFIED 


SYSTEM  TRANSFER  FUNCTION 

POLES  REAL 

1  -1.0 

2  -10.0 

3  -20.0 

4  -33.0 

5  -100.0 


IMAGINARY 
0.0 
0.0 
0.0 
0.0 
0.0 


J 
J 
J 
J 
J 


ZEROES 

1 


REAL 
-5.0 


IMAGINARY 
0.0 


GAIN  CONSTANT  =  2640000.0 


TESTS  FOR  EXAMPLE  10 


RUN 

NO. 

OF    POLES 

a 

5 

b 

3 

c 

2 

NO.  OF  ZEROES      RESULT 

1  See  Table  of  Data 

0  See  Table  of  Data 

0  See  Table  of  Data 


INPUT  FUNCTION  SIMILAR  TO  THE  ONE  SHOWN  IN  FIG.  5. 
INITIAL  CONDITIONS  EQUAL  ZERO. 

SEE  FIG.  6  COMPARING  THE  RESPONSES  OF  THEORETICAL 
AND  IDENTIFIED  SYSTEMS  TO  A  UNIT  STEP  INPUT. 
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IDENTIFICATION  OF  UNKNOV7N  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

•   1 

-0.94505946 

0.0 

J 

2 

-10.807637 

0.0 

J 

3 

-32.034223 

0.0 

J 

4 

-62.693979 

0.0 

J 

5 

-120.97375 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 

-4.9340575 

0.0 

J 

GAIN  CONSTANT  =  2741280.0 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

248.538.7 

2 

2992971.5 

3 

403446.21 

4 

16013.272 

5 

227.45465 

B  VECTOR 

5 

2741280.0 

C  VECTOR 

1 

4.9340575 

2 

1.0000000 

100 


REMARKS: 

ASK  FOR  FIVE  POLES  ONE  ZERO. 

RESULTS  NOT  VERY  GOOD  FOR  THE  HIGHER  FREQUENCY  POLES. 

SEE  FIG.  6  COMPARING  THE  RESPONSES  OF  THEORETICAL  AND 
IDENTIFIED  SYSTEMS  TO  A  UNIT  STEP  INPUT. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (b) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -2.6408186 

2  -24.059549 

3  -24.059549 


ZEROES 


REAL 


IMAGINARY 

0.0 

J 

0.0 

J 

0.0 

J 

IMAGINARY 

GAIN  CONSTANT  =  10410.082 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 
2 
3 

B  VECTOR 

3 
C  VECTOR 

1 


3331.3020 
1388.5395 
50.759916 

10410.082 

1.0000000 


REMARKS : 

ASK  FOR  THREE  POLES  NO  ZEROES. 

RESULTS  CAN  BE  CONSIDERED  A  LOWER  ORDER  MATHEMATICAL 
MODEL . 

SEE  FIG.  6  COMPARING  THE  RESPONSES  OF  THEORETICAL  AND 
IDENTIFIED  SYSTEMS  TO  A  UNIT  STEP  INPUT. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (c) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -5.6608772 

2  -13.680660 


IMAGINARY 

0.0  J 

0.0  J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  182.89790 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 

1  77.444535 

2  19.341537 
B  VECTOR 

2  182.89790 

C  VECTOR 

1  1.0000000 

REMARKS : 

ASK  FOR  TWO  POLES  NO  ZEROES. 

RESULTS  CAN  BE  CONSIDERED  A  LOWER  ORDER  MATHEMATICAL 
MODEL. 

SEE  FIG.  6  COMPARING  RESPONSES  OF  THEORETICAL  AND 
IDENTIFIED  SYSTEMS  TO  A  UNIT  STEP  INPUT. 
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12  3  4 

TIME-SECONDS 

UNIT  STEP  RESPONSE  OF: 

(a)  System  to  be  identified 

(b)  Fifth  order  identified  model 

(c)  Third  order  identified  model 

(d)  Second  order  identified  model 

FIGURE  6 
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EXAMPLE  11 


SYSTEM  TO  BE  IDENTIFIED 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -1.0 

2  -10.0 

3  -33.0 

4  -100.0 


ZEROES 
1 


REAL 
-5.0 


IMAGINARY 

0.0 

J 

0.0 

J 

0.0 

J 

0.0 

J 

IMAGINARY 

0.0 

J 

GAIN  CONSTANT  =  52500.0 


TESTS  FOR  EXAMPLE  11 


RUN 

NO. 

OF  POLES 

a 

4 

b 

3 

c 

2 

NO.  OF  ZEROES      RESULTS 
1  Good 

0        See  Table  of  Data 
0        See  Table  of  Data 


INPUT  FUNCTION  IS  SIMILAR  TO  THE  ONE  SHOWN  IN  FIG.  5. 
INITIAL  CONDITIONS  EQUAL  ZERO. 

SEE  FIG.  7  &  8  COMPARING  THE  RESPONSES  OF  THEORETICAL 
AND  IDENTIFYING  SYSTEMS  TO  A  STEP  INPUT. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (a) 


SYSTEM  TRANSFER  FUNCTION 


POLES 

REAL 

IMAGINARY 

1 

-1.0182594 

0.0 

J 

2 

-10.853145 

0.0 

J 

3 

-35.841486 

0.0 

J 

4 

-108.22539 

0.0 

J 

ZEROES 

REAL 

IMAGINARY 

1 
GAIN 

CONSTANT 

-5.0491588 
=  64711.781 

0.0 

J 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 

42867.599 

2 

47640.814 

3 

5600.2860 

4 

155.93828 

B  VECTOR 

4 

64711.781 

C  VECTOR 

1 

5.0491588 

2 

1.0000000 

REMARKS: 

ASK  FOR  FOUR 

POLES  ONE  ZERO. 

RESULTS  GOOD. 
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SEE  FIG.  7  &  8  COMPARING  THE  RESPONSES  OF  THEORETICAL 
AND  IDENTIFIED  SYSTEMS  TO  A  UNIT  STEP  INPUT. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (b) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -31.063053 

2  -31.063053 

3  -2.9265019 


IMAGINARY 
-36.791868 
36.791868 
0.0 


J 
J 
J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  2  8337.477 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 


1 
2 
3 

B  VECTOR 

3 
C  VECTOR 

1 


6785.2550 
2500.3670 
65.052607 

28337.477 

1.0000000 


REMARKS : 

ASK  FOR  THREE  POLES  NO  ZEROES. 

RESULTS  CAN  BE  CONSIDERED  A  LOWER  ORDER  MATHEMATICAL 
MODEL . 

SEE  FIG.  7  &  8  COMPARING  THE  RESPONSES  OF  THEORETICAL 
AND  IDENTIFIED  SYSTEMS  TO  A  UNIT  STEP  INPUT. 
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IDENTIFICATION  OF  UNKNOWN  SYSTEM  -  Run  (c) 


SYSTEM  TRANSFER  FUNCTION 
POLES  REAL 

1  -3.4314173 

2  -36.703534 


IMAGINARY 

0.0  J 

0.0  J 


ZEROES 


REAL 


IMAGINARY 


GAIN  CONSTANT  =  537.30396 

SYSTEM  STATE  VARIABLES  (PHASE  FORM) 
A  VECTOR 

1  125.94514 

2  40.134951 
B  VECTOR 

2  537.30396 

C  VECTOR 

1  1.0000000 

REMARKS : 

ASK  FOR  TWO  POLES  NO  ZEROES. 

RESULTS  CAN  BE  CONSIDERED  A  LOWER  ORDER  MATHEMATICAL 
MODEL . 

SEE  FIGs.  7  &  8  COMPARING  THE  RESPONSES  OF  THEORETICAL 
AND  IDENTIFIED  SYSTEMS  TO  A  UNIT  STEP  INPUT. 
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UNIT  STEP  IRESPONSE  OF: 

|  (a)  System  to  be  identified 

j (b)  Fourth  order  identified  model 

I (c)  Third  order  identified  model 

i  (d)  Second  order  identified  model 


2  4 

TIME-SECONDS 

FIGURE  7 
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UNIT  STEP  RESPONSE  OF : 
(a)  System  to  be  identified 

Fourth  order  identified  model 


(b) 
(c) 
(d) 


Third 


order  idejntified  model 


Second  order  identified  model 


2  4 

TIME-SECONDS 

FIGURE  8 
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IV.   CONCLUSIONS 

The  method  of  multiple  integrations  is  a  practical  method 
for  identifying  lumped  parameter,  linear,  time  invariant 
systems.   It  works  well  in  the  absence  of  noise  up  to  a  fourth 
order  system  if  the  input  function  is  carefully  selected. 
For  higher  order  systems  the  identification  is  in  general  not 
good,  but  an  approximation  of  the  system  can  be  obtained. 
If  desired,  the  method  can  also  be  used  to  obtain  a  reduced 
order  mathematical  model  of  the  system,  after  the 
identification. 

The  accuracy  of  an  identification  depends  very  strongly 
on  the  noise  present.   Small  R.M.S.  values  of  noise  are 
enough  to  make  the  identification  bad.   Noise  with  a  R.M.S. 
value  of  about  0.5  volts  and  below  do  not  practically  affect 
the  identification,  when  inputs  like  the  ones  ued  in  this 
research  are  present,  independently  of  where  the  noise  is 
present  -  input,  output  or  both  -.   Accuracy  also  depends 
very  much  on  the  type  of  input  function  used  to  drive  the 
system. 

The  computational  requirements  of  the  method  are  not 
excessive,  however,  an  accurate  A/D  converter  system  is  very- 
desirable  and  the  better  this  operation  is  done,  the  more 
accurate  will  the  identification  be,  since  the  accuracy  of 
an  identification  is  usually  comparable  to  the  number  of 
exact  significant  digits  in  the  input-output  data. 
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The  limitation  of  the  method  to  low  order  systems  is 
primarily  due  to  the  algorithm  used  to  solve  the  over-deter- 
mined set  of  linear  equations.   As  better  algorithms  become 
available,  it  will  be  possible  to  identify  higher  order 
systems. 

The  results  obtained  when  reduction  of  the  system  order 
was  desired  were  not  very  good  at  first  sight,  but  when  the 
steady  state  gain  of  the  lower  order  models  was  adjusted, 
making  it  equal  to  the  one  of  the  non-reduced  identified 
model,  the  approximation  can  be  considered  good.   However, 
when  a  pair  of  dominant  complex  roots  is  present  in  the 
system  a  better  reduced  order  model  should  be  expected. 
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APPENDIX  I 

A.   OBTAINING  THE  INPUT  -  OUTPUT  RECORD 

The  identification  technique  described  in  Chapter  II  is 
based  on  knowledge  of  the  input  and  output  record  and  since 
the  IBM  360  Digital  Computer  was  used  to  perform  the  inte- 
grations of  equation  (5) ,  the  record  should  be  available  in 
a  digital  form.   Thus,  and  A/D  (analog  to  digital)  operation 
should  be  performed.   This  A/D  operation  was  accomplished 
using  the  Hybrid  Computer  System  (XDS  9300/Ci  5000)  located 
in  the  Electrical  Engineering  Computer  Laboratory.   The 
output  of  this  computer  is  a  seven  track  magnetic  tape  in 
which  the  input  and  output  functions  are  recorded  in  digital 
form.   Since  the  XDS  9300/Ci  5000  system  is  used  on  a  self- 
service  basis,  a  brief  discussion  of  what  was  done  to  obtain 
the  record,  is  given  here,  so  that  the  same  procedure  could 
be  used  in  the  future  whenever  an  A/D  conversion  is  needed  - 
it  should  be  pointed  out  here  than  the  computer  center  staff 
was  always  predisposed  to  assist  in  the  operation  of  the 
system. 

The  analog  part  of  the  system,  the  Ci  5000,  has  two  patch- 
boards which  are  the  heart  of  the  system.   The  logic  patch- 
board was  utilized  to  set  the  sampling  rate,  at  which  the 
analog  signal  was  sampled  -  10KHZ  in  this  particular  case. 

This  sampling  rate  is  the  maximum  that  can  be  achieved 
with  this  computer,  and  can  be  changed  (decreased)  by 
altering  a  counter-like  number  display  located  in  the 
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frontal  panel  of  the  computer  -  FPOO  and  FP04.   Also  in  this 
patchboard  the  base  frequency  -  1KHZ  -  used  to  determine  the 
moments  at  which  the  input  function  is  changed  step-wise,  is 
set.   The  same  counter-like  number  display  can  be  used  to 
change  this  frequency  -  FF10  and  14.   In  the  digital  patch- 
board, the  wiring  that  controls  the  digital  part  of  the  A/D 
operation  is  also  set. 

The  analog  patchboard  was  used  to  introduce  the  proper 
isolation  between  the  different  components  of  the  different 
systems  tested  and  to  provide  the  desired  amplification 
factor.   This  patchboard  was  also  used  in  the  generation  of 
the  input  signal.   This  signal  was  obtained  in  the  following 
way.   An  oscillator  generates  a  sine  wave  in  which  amplitude 
and  frequency  can  be  adjusted  by  two  potentiometers  -  P043 
and  P417  -;  the  .sine-wave  is  amplified  and  passed  through  a 
zero  order  hold  circuit  which  output  is  a  step-wise  wave  form, 
that  was  used  as  the  input  function.   The  sine  wave  was  also 
used  to  control  the  beginning  of  the  sampling.   The  wiring  of 
both  patchboards  was  set  by  the  computer  center  staff  and 
those  patchboards  -  numbered  21  and  21  -  are  reserved  for 
the  A/D  operation. 

The  input  and  output  signals  are  fed  from  the  analog 
part  of  the  computer  system  to  the  digital  part.   The 
Digital  part  interfaced  with  the  analog  is  a  Xerox  Data 
System  -  XDS9300. 

Two  tapes  drive  units  are  available  as  input-output 
devices.   One  printer  is  used  for  program  listings  and 
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21.  With  the  sense  switches  at  the  XDS9300  console  off, 
depress  CLEAR  and  CLEAR  FLAGS  buttons  simultaneously. 

22.  Mount  the  magnetic  tape  following  the  instructions 
available  on  the  inside  of  the  protective  doors. 

23.  Select  density  556  bits  per  inch  with  the  correspond- 
ing knob  on  the  face  of  the  tape  driver. 

24.  Read  in  the  A/D  program,  described  in  Chapter  III 

and  listed  in  Chapter  IV,  following  the  instructions  available 
on  the  console  of  the  XDS9300. 

25.  Enter,  through  the  teletype,  the  parameters  to  be 
used  for  the  A/D  program.   This  program  has  seven  options  or 
instructions  which  are  entered  also  through  the  teletype  by 
a  single  digit  number  and  depressing  the  carriage  return 
button.   If  a  further  response  is  required  the  computer  would 
indicate  the  format.   The  options  are: 

1  =  Desire  to  enter  new  parameters. 

2  =   Start  digitizing  the  analog  input  signals 
(normally  actuated  by  manual  switch  operation  on  Ci  5000 
control  panel;  (DS1  up  and  key  COMPUTE  depressed) . 

3  =   Write  an  END  OF  FILE  on  the  magnetic  tape. 

4  =   Rewind  the  tape  to  the  load  point  (beginning) . 

5  =   Skip  the  following  number  of  END  OF  FILE  marks. 

6  =   Print  the  following  amount  of  digitized  data. 

7  =  Actuate  the  A/D  program  for  the  next  single 
block  of  data  encountered  and  feed  the  signal  to  the  strip 
chart  recorder  (which  must  be  running  prior  to  this  option 
selection  connected  to  the  analog  patchboard  on  the  Ci  5C00. 
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The  rest  of  the  counters  should  be  set  to  zero. 

7.  Depress  the  master  clock  control,  RUN. 

8.  Set  DFOO  (Delay  flip-flop)  to  .1  ms.  and  adjust  it, 
until  its  corresponding  light  on  the  front  panel  of  the 

Ci  5000  comes  on. 

9.  Set  DF06  in  .1  ms.,  DF07  in  1  sec.  and  DF10  in  1  ms. 

10.  Set  all  the  DS  (digital  switches)  to  the  center 
position. 

11.  Turn  on  the  TRC. 

12.  Set  switches  XI  to  the  left  and  X2  to  the  right. 

13.  Select  amplifier  042  in  the  Ci  5000  console,  adjust- 
ing P417,  DF06  and  F10,  F14  to  obtain  the  desired  input 
function  which  can  be  observed  in  the  TRC  when  the  DS1 
switch  is  set  in  the  up  position  and  the  button  COMPUTE 

is  depressed  in  the  Ci  5000  console. 

14.  Set  the  R.T.  clock  switch  on  XDS9300  to  EXIT. 

15.  If  the  XDS9300  is  not  energized,  depress  simultaneously 
RESET  and  POWER  buttons. 

16.  If  the  teletype  is  not  ON,  energize  it  by  turning  the 
knurled  knob  to  ON. 

17.  If  hte  Line  Printer  is  not  ON,  depress  the  POWER 
switch  -  indicator  on  the  face  of  the  printer.   About  one 
minute  later  the  lower  half  of  the  Power  button/lights  up 
indicating  the  printer  is  energized. 

18.  Depress  READY  SWITCH  INDICATOR. 

19.  Depress  the  POWER  ON  button  on  the  card  reader. 

20.  Depress  IDLE  button  on  the  console  of  the  XDS9300. 
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outputs.   A  card  reader  is  associated  with  the  digital  com- 
puter as  the  avenue  of  input  for  compiling  programs.   A 
teletype  unit  is  used  for  parameter . inputs  and  instructions 
to  the  digital  computer  program.   A  control  console  displays 
information  of  the  computer  and  has  the  corresponding  control 
switches.   A  recorder  is  also  available  with  which  the 
records  of  the  input  and  output,  after  the  A/D  operation  is 
finished,  can  be  recorded  in  a  strip  chart. 

After  the  preceding  general  description  of  the  system, 
the  following  procedure  can  be  followed  to  take  a  record. 

1.  Set  the  analog  and  digital  patchboards  in  the 
Ci  5000. 

2.  Connect  the  output  of  the  amplifier  042  to  the 
input  of  the  system. 

3.  Connect  the  output  of  the  system  to  the  input  of 
the  amplifier  036. 

4.  If  the  Ci  5000  is  not  energized,  depress  on  switch 
indicator  at  lower  left  corner  of  display  panel  of  Ci  5000. 
Depress  KEY  BOARD  and  POTSET  buttons  on  the  display  panel. 

5.  Enter  the  maximum  value  of  the  sine  wave  through 
the  POTSET,  POT  and  SERVO  keys. 

6.  Set  the  counters  FF ' s  in  the  following  positions: 
FF00  in  0,  FF04  in  9,  for  10KHZ  sampling  rate. 
FF10  in  2,  F14  in  4,  for  40HZ  rate  of  changing, 

step-wise,  the  input  function. 

FF30  in  0,  if  only  one  record  per  run  is  desired; 
increasing  FF30  the  number  of  records/run  can  be  adjusted. 
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The  parameters  to  be  entered  in  point  25  are: 

NSAMP   =   Total  number  of  samples  to  be  taken  from 
each  channel  -  4002  in  this  case. 

NCHAN   =   Number  of  channels  to  be  digitized  -  2  in 
this  case. 

NREC   =  Number  of  records  to  be  digitized  per  run; 
this  number  should  be  equal  to  the  one  set  in  the  counter  FF30 
(Ci  5000)  plus  1;  in  this  particular  problem  NREC  =  1. 

ITAPE  =   The  identifying  number  of  the  tape  drive  in 
use.   This  number  was  set  when  the  tape  was  mounted. 

NDEL   =   Controls  the  delay  between  the  digital  and 
analog  signal  in  option  7 . 

After  entering  the  parameters  the  *  character  should  be 
typed  on  the  teletype  followed  by  RETURN.   At  this  point 
everything  is  ready  to  take  a  record;  by  typing,  2  -  RETURN, 
on  the  teletype,  and  setting  switch  DS1  -  UP  and  pressing 
COMPUTE  on  the  Ci  5000.   After  the  record  was  taken,  write 
-  3  -  (END  OF  FILE) .   If  it  is  desired,  the  record  can  be 
recorded  in  a  strip  chart  recorder. 

B.   SEVEN  TO  NINE  TRACK  TAPE  CONVERSION 

In  the  magnetic  tape  generated  in  point  A,  the  input- 
output  record  is  written  in  a  seven  track  tape,  but  since 
this  tape  is  going  to  be  used  with  the  IBM  360  Digital 
Computer  and  this  computer  uses  NINE  TRACK  TAPES,  a 
conversion  of  the  seven  track  tape  into  a  nine  track  tape 
is  required. 
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This  operation  is  done  by  the  IBM  360,  using  the  program 
CONVERT,  which  is  described  in  the  Technical  Note  No.  0211-05 
second  edition  -  PROCEDURES  FOR  CONVERTING  7 -TRACK  MAGNETIC 
TAPES  TO  9-TRACK  MAGNETIC  TAPES  by  Sharon  D.  Raney  [21] . 
This  description  is  available  at  the  Computer  Center  in 
Ingersol  Hall.   A  listing  of  the  CONVERT  program  is  given 
in  Chapter  IV. 
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Several  practical  examples  are  presented  to  demonstrate  the 
accuracy  and  present  capabilities  of  the  procedure. 
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